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Abstract 
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1 INTRODUCTION 


It is generally accepted that QCD generates two phase transitions, a decon- 
hnement one, corresponding to a transition of a hadron gas to a quark-gluon 
plasma and a chiral one, corresponding to the passage from a chiral symmetry 
breaking phase to a phase where chiral symmetry is restored. Lattice QCD 
calculations, based on the estimate of the Polyakov loops, show that the two 
transitions coincide in the pure gauge case. Lattice calculations including mas¬ 
sive quarks [1] (albeit not yet in a perfect shape) seem to indicate that this 
situation persists [2-4] after the introduction of quarks. However the order of 
the transition is changing, from hrst order, for the pure gauge case, to second 
order with reasonable values of the quark masses. The case of the usually ac¬ 
cepted values for the bare u, d and s masses seems to be close to the boundary 
separating first and second order transitions at zero quark chemical potential. 
Nothing is really known for non-zero chemical potential. 

Whereas the deconfinement cannot be studied in existing satisfying models, 
the features of spontaneous chiral symmetry breaking and of its restoration 
can be studied in the cr-model or in the Nambu Jona-Lasinio (NJL) model [5]. 
The latter seems to be very successful in the description of many features of 
hadronic physics in the non-perturbative domain of QCD, despite of the ab¬ 
sence of confinement and of renormalizability. It seems that the chiral aspects 
are much more important than the confinement features, for static properties 
at least. This model is thus well suited to study the chiral transition of an in¬ 
finite uniform quark system, and the thermodynamics of the phases. The pure 
NJL model with SU(2) flavour has been used for this purpose by many authors 
[6-11]. Depending upon the unavoidable choice for free parameters left in this 
kind of model, the transition can be of first or second order [12], with a critical 
temperature of the order of 150-200 MeV at zero chemical potential. Most of 
these approaches are restricted to the classical (mean field) approximation. 
Let us nonetheless notice the works of Zhuang et al. [11] and of Blascke et ah 
[13], who calculated the thermodynamic potential to the order where W 
is the number of colors, which amounts roughly to include the meson degrees 
of freedom. The latter seem to play an important role at low temperature 
only and does not affect significantly the position and the properties (order, 
discontinuities,...) of the transition [11]. The chiral transition has been studied 
for the SU(3) flavour case at zero chemical potential [14-17] and at non-zero 
light quark {u,d) chemical potential [18]. 

Here, we want to make a systematic study of the chiral phase transition fea¬ 
tures for the SU(3) flavour case for any value of the light quark (taken 
equal to p^) and of the strange quark chemical potentials. For this, we 
use an extension of the NJL model, which implements the trace anomaly of 
QCD through the introduction of a scalar dilaton field, and the axial anomaly. 
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The former allows the description of the gluon condensate and the latter is 
dictated by the low energy hadron phenomenology. The model as well as its 
predictions for the temperature and density dependences of the meson masses 
are described in ref. [17]. We will evaluate the thermodynamic potential in the 
mean field approximation within this model and make a systematic survey of 
the transition in the T, and p* spaces. The main thermodynamic quantities 
will be derived for each phase. We use the assumption of a vanishing current 
mass for the light quarks. The non-zero current s-quark mass constitutes the 
only flavour SU(3) symmetry breaking. 

The paper is organized as follows. The NJL model used here is schematically 
described in section 2, as well as the associated action. Section 3 is devoted 
to the localisation of the chiral phase transition in the T, /r„ and /r* space 
for two choices of the free parameters of the model. For the first one, the 
transition is only of first order, whereas the second one allows the possibility 
of a second order transition. Section 4 deals with the pressure, energy and 
entropy densities in the different phases. The limiting values in the chiral 
phase are particularly examined. Section 5 contains our conclusion. 


2 THE A-SCALING NJL MODEL 


The model used in this paper is extensively described in Ref. [17] under the 
name “A-scaling model”. We just recall some of the tools useful to the under¬ 
standing of the present paper. 

We start from the SU(3) effective Euclidean action : 


hff (v?, X, yw, /?) = In 

+ J + j d‘^xC^ + j d'^xCu^ii) (1) 

where the quark fields have been integrated out. The meson fields write : 

T’a = (c^a,7ra) , F^ = (A^, WgAa) a = 0, ...,8 (2) 

where A^ are the usual Gell-Mann matrices (^Aq = \/|l) the gluon con- 

densate is represented by the scalar dilaton field y oc . The lagrangian 

given by 
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takes account of the QCD trace anomaly via the term [22,23]. The 

Lhi(l) axial anomaly is introdnced via the mass term : 


^Ua{1) 



( 4 ) 


In Eq. (1), the trace has been regularized by an ultra-violet cut-off which is 
X-dependent in order to preserve the scale invariance of the quark loop term. 
The chemical potentials are introduced via the vector field W\ = 0^, 

where /r is a shorthand notation for fi = diag fig), if one neglects the 

isospin breaking (/i^ = fid). In the same way the quantity m represents the 
diagonal matrix m = diag {mu, mu, mg) in flavor space. 


We work in the classical approximation so that the partition function writes : 




exp 


^eff {^u, Mg, Xs 


( 5 ) 


where (Mg) and Xs denote respectively the constituent u-quark (s-quark) 
mass and the glnon condensate in the hot and dense system and are thus 
[I- and /^-dependent {{3 = 1/T). In the vacuum, the corresponding quantities 
are noted M°, and xo- These three quantities (in the vacnnm or in the 
system) minimize the thermodynamical potential or equivalently the action 
(1), yielding the value entering eq. (5). The relationship between Mu, Mg 
and the stationary fields cTq, (j| is given in [17]. We write : 


4°// (A4, X.) = (M„, M„ X,) + (M„, M„ x. 


(S) 


where 

(^n, Mg, Xs) =-Tr In {-i^ - l;T-f m + ra(pa)‘' 

+ Tr In {—i^m-\-Ta(pa) , (7) 

where the index c indicates that the traces have to be evaluated with the 
stationary values of the fields. The contribution (7) to the action is finite and 
writes : 


2N r 

1 ( 1 , 0 } (Mu, M„ X.) = - '£<^.1 1 = 

M j—u.S 


'‘dk 


X {In (1 -F exp [-j3 {Ei + fii)]) -|- In (1 -|- exp [-/? (Ej - /Tj)])} 


( 8 ) 


with Ei = y k‘^ + = 2, = 1 and Q the volume of the system. The 

second term of the r.h.s. of (6) corresponds to fiu = fig = 0, (3 = ^ = oo. It 
writes : 
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+ 



^(0,oo) (Mu, Xs) = 




-TtAx, + m + VaifaT 




( 9 ) 


which yields with a d'^k cut-off : 


2Nr 




X 


{[(Ax.) 

-\-PQ 


-Mi 


In 


(^Xs)^ + Mi 


l=U,S 
2 i\ ^ 2 


-{AxsY-m, 


{AXsf + MY In Mi 




E o.iM^H^ 


mi\‘^ 

m) +16 


x/ln 


,Xg, 


{xt - Xa) 


( 10 ) 


Our model contains seven parameters : the bare quark masses (m„,ms), the 
strengths (a^,6^), the vacuum gluon condensate xo, the cut-off A and the 
parameter ^ of the axial anomaly. Five of these seven parameters are fixed to 
reproduce the vacuum masses of the pion tUtt = 139 MeV, of the kaon m-K 
= 496 MeV, of the rj' m,,/ = 958 MeV and of the glueball mcL = 1-3 GeV, 
together with the pion decay constant /n- = 93 MeV. The two other parameters 
remain free and they have been chosen to be the vacuum constituent u-qnark 
mass (related to a^) and the vacnum gluon condensate xo- We give results 
for two sets of parameters : 

M° = 600 MeV , xo = 125 MeV (11) 

M° = 300 MeV , Xo = 80 MeV. (12) 


This choice is driven by the fact that at /i„ = = 0, the phase transition is 

of first order (Tc = 168) MeV with the first set and of second order (Tc=150 
MeV) with the second set. Note that in the SU(2) case [12], Tc was fonnd to 
be 140 MeV with the second set, in agreement with the lattice data [19]. We 
then investigate how finite density can affect the order of the transition. Since 
we are interested in the determination of the parameters (Tc, /r„c, A^sc) at the 
phase transition, we choose to work in the limit = 0. 
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3 CHIRAL PHASE TRANSITIONS 


The quark condensates are related to the constituent masses in the following 
way : 


^ (A4, A/,.x.) = -i<.q;Af„ 




(13) 

(14) 


The coupling between these two condensates comes through their coupling to 
the gluon condensate The value of is indeed coupled to those of and 
Mg by the search of the minimum of the action ( 6 ) (see Eq. ( 1 )). This coupling 
can be better visualized by looking at the “three coupled gap equations” (see 
[17]) that amount to minimize the action. 

The chiral phase transitions happen when the quark condensates go to zero. 
There is a priori no reason why the two condensates should vanish at the same 
time. This is illustrated in Fig. 1 where the open circles represent the phase 
transition jigc versus jiuc, corresponding to a vanishing u-quark condensate at 
Tc = 0 (a), Tc = 100 MeV (b), Tc = 161 MeV (c), for the set of parame¬ 
ters ( 11 ). The transition for the s-quark is given by the crosses. At sulRciently 
small temperatures, one observes that the chiral transition of the u- and s- 
condensates are independent of one another (as in a pure NJL model without 
’t Hooft determinant), except in a small region of the space (/iso hue) where 
the coupling via the gluon condensate is sufficiently strong to force the two 
quark condensates to go to zero at the same time. When the temperature is 
high enough, only this region of strong coupling is present. Whatever the crit¬ 
ical point of the three-dimension space {T^ Huc, fJ^sc), the transition is of first 
order. This is iliustrated in Fig. 2 in the case of strong coupfing (between the 
condensates) (a) (Tc ~ 120 MeV, /i„c = /^sc = 250 MeV) and of weak coupfing 
(b) {T^c = 140 MeV, Tgc = 160 MeV, = 200 MeV, fisc = 0 MeV). Within 
the present modef, a discontinuity (transition) of a quark condensate is always 
accompanied by a discontinuity of the gluon condensate, in contradistinction 
with the suggestion of Ref. [20]. Note that y 7 ^ 0 above the transition. Fig. 3 
summarizes the previous results in the form of an artistic view. The space is di¬ 
vided into four regions : the inner one corresponds to quarks in their condense 
phase ; the outer one can be seen as the quark-gluon plasma, the two remain¬ 
ing regions (extending along the axes to infinity) correspond to mixed phases 
with free u (s)-quarks together with s (u)-quarks in the hadronic phase if one 
assumes that the chiral transition corresponds to deconfinement. A remarkable 
prediction of the present model has to be singled out : if non-strange matter 
at sufficiently large up chemical potential is heated, the matter will transform 
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first to a mixture of free light quarks surrounding strange hadrons (with total 
vanishing strangeness) before being transformed into a quark gas. 


The densities and Ps are related to the chemical potentials and temperature 
in the following way 


Pi 



dllfA 

dPi I 


7J.2 


oo 

j (n+j — n_j) dk i = u,s (15) 


where 


n±i = - 7 — / -^ i = u,s. (16) 

1 + exp |/? ± pij I 


In Eq. (15), takes into account the d-quark contribution via a„ = 2. When 
the transition is of first order, there is a jump in the density at the transition. 
Just before the transition, the p^’s are given by (15) and (16), with Mf ^ mf 
{Pi bef), while just after the transition, the densities (p, aft) increase up to 
values corresponding to Mi = mi. A drawback of the present model is that at 
Tc = 0, pi bef = 0, so that the system can not exist in the hadronic phase. As 
Tc departs from zero, pi bef remains small. As an example, for puc = Psc = 250 
MeV, one has = 0.362, = 0.050 and = 2.646, = 0.978 

where po = 0.51 fm~^. Within the present model, it is practically impossible 
to eliminate this drawback, as the only way would correspond to having small 
constituent masses in the hadronic phase while keeping large gluon condensate. 
Such a choice of parameters is not suitable [12] since it would yield too high 
Tf.. A hope lies however in the inclusion of vector mesons which are expected 
to shift the critical densities up to acceptable values without modifying the 
qualitative features of the phase transition [25], as the latter are basically 
determined by the chiral mesons and their coupling to the condensates. 

The previous results are not qualitatively changed drastically when shifting 
to the other choice of the parameters (Eq. (12)), except for the order of the 
transition which can become of second order at sufficiently low chemical po¬ 
tentials. This is illustrated in Eig. 4a. At Pu = Ps = 0, the u-quark condensate 
vanishes at T^c ~ 150 MeV (140 MeV in SU(2)) while chiral transition hap¬ 
pens at Tsc ~ 200 MeV for the s-quark. At higher chemical potentials, the 
transition becomes of first order but the chiral transitions for the u- and s- 
quarks remain decoupled (see Fig. 4b). Note, however, that a change of the 
gluon condensate is always associated with the transition of any quark con¬ 
densate. Figs. 5 and 6 are the analogs of Figs. 1 and 3, respectively for the 
set (12). Figs. 1 and 5 are qualitatively similar except at high temperature 
where the two quark condensates are coupled whatever the couple {puc, Psc) 
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for the set (11) while the ti-quark can only exist in the “deconfined” phase 
with a chiral transition at ~ 220 MeV whatever Huc for the strange quark, 
in the case of set (12). The results presented in Fig. 5 can be visualized in the 
3-dimension space by use of the artistic view of Fig. 6. 


4 THERMODYNAMICS 


We now turn to the study of the thermodynamics of a hot and/or dense system. 
We will mainly focus on the behaviour with respect to the temperature of the 
pressure, the energy density and the entropy density. These three quantities 
are respectively given in the classical approximation by : 


P = 


iain(Ze/Z0) 




dll^^ {Mu,Ms,Xs) a/^0,0,) (M°,M°,xo) 






fs dn 

^ (M„, A/„ X.) + (M„. A/„ x.) - Afo.oo) (A/„", A/.°, X'o) 


fXt L 


(17) 





+ 


pn L 


%,«,) (A4, M., x.) - (Af.“, Af», xo); 


(A4“. Af,“, xo) 


(18) 


_S „ E „ 

' l=U,S 


0 




(19) 


where we have substracted the vacuum contribution from the partition 
function in order to ensure vanishing pressure, energy and entropy in the 


vacuum. 



Following the same prescription as in [21], we can write 

P Pideal gas P l^U) 

^ ^ideal gas P {P j /^s) 

Pideal gas P ^ideal gas ^ ^ l-^iPi 


( 20 ) 

( 21 ) 

( 22 ) 


with 


Pi 


ideal gas 


11,0} (M... M., X.) _ N, ^ ,23) 


fsn 


Stt^ 


l = U,S Q 


\ ( Q 1 _^ ^ 

^ideal gas = -^ ^ ^ ~ ^ X/ j "^S’ ^s) 


Nc f 

^ ^ J k Ei ij^+i P 'IT'—i) 


dk 


(24) 


I,S Q 


and B denotes the (temperature and density-dependent) bag pressure 


B W. f,u, ft) = P {/&,„) (M„, M„ X.) - /5,„) (Ml Ml Xo) } . (25) 


(3^ 


The definition (25) differs from the one used in [17,22,23] where the bag pres¬ 
sure is defined as: 

B = P {y„,„, (0,0.0) - (Ml, Ml Xo) } 

= ^”icLXl- (26) 

The temperature behaviour of P is exhibited in Fig. 7a, for the two sets of 
parameters (11) and (12) in the case /i„ = /i^ = 0. Due to the non-vanishing 
value of the strange quark bare mass mg, the temperature behaviour above Tc 
is not exactly the usual law in T^. One can write (see Appendix A) 


Pideal gas ' 


7 N 

‘ AT 2 t -.4 -'’c 2 t-t 2 

— NcTT T - mzT — 

60 12 * 


Nc 4 1 

s?’"* 


s 

~nf 


+ 


Nc 


-27) 




(27) 
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where 7 is the Euler constant. At the same time, B becomes T-independent 


B 


1 



(0,ms,Xs) 


-I 


C 

(0,00) 


K“.ACxo)}, 


( 28 ) 


where is the stationary value for T > Tc. Fig. 7a shows an excellent 
linear fit of P above chiral restoration : 

P ^ 3.40 - 2.097 10■^ (29) 


for the set of parameters ( 11 ) and 

P ^ 3.40 - 1.195 10■^ 


(30) 


for the set (12). The coefficient of is not exactly equal to the coefficient of 
the same power in (27), which is equal to 3.454. In fact, due to the limited 
range of investigated in Fig. 7, the remaining terms in Eq. (27) change 
the apparent slope and the constant term of the curve. When a fit with the 
explicitly given terms of Eq. (27) plus an adjusted constant C is performed, 
the latter turns out to be ~ 207 MeV for set (11) and (7^/^ ~ 179 
MeV for set (12). These values are close to 209 MeV (set (11)) and 

~ 183 MeV (set (12)) obtained from a direct calculation using Eq. (28). 
This underlines the consistency and the precision of our numerical results 
as well as the rapid convergence of expansion (27), although the expansion 
parameter rng/T is not very small. Moreover, note that the numerical values 
of 5^/^ given by (28) or extracted from the fits are not so far from the values 
of B^/'^ given by (26) (S^/^ = 201 MeV (set 11), 161 MeV (set 12)). 


Below the phase transition, the temperature behaviour depends on the order 
of the transition. With the set (11) (first order), (3Mi {i = u,s) remain large, 
so that the T behaviour is governed by the law 


P 


ideal gas 


(27Tf^ 




(31) 


It should be noted that the masses Mj are decreasing with temperature, which 
makes the variation of the pressure more rapid that the factor in Eq. (31). 
The masses Mj being larger with set (11) than with set (12), the pressure 
increases more slowly in the first case and the change of slope at the transition 
is more pronounced. 

The law for Pideai gas does not contain the explicit gluon degrees of freedom 
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which would yield : 


2 'v 

Pid^ai gas ~ ^ [^N g + -N,Nf4^ ^ 5.2 T\ (32) 

A first step in the inclusion of the gluon thermodynamics consists in going 
beyond the classical approximation by adding one loop corrections. One then 
has to add 


tt) 

1 I" dk 1 

+ (^1 — exp/^-^fc^ + my 

where the sum runs over the 19 mesons : (9+1) scalar and 9 pseudo-scalar 
mesons. When chiral symmetry is restored, all the meson masses {rrij) (except 
for the dilaton) increase to infinity for increasing temperature so that the 
product (drrij remains large, and : 



P, 


one loop 




(2 


In 


1 — exp + 


rud 


Pone loopimeson) 




dk 


k^ 


27r" 


^-5/2 




( 27 r) 


3/2 




ZI2 -Pmj 


(34) 


in the non-relativistic limit. The exponential factor makes the one loop correc¬ 
tion terms due to the 18 mesons vanishingly small compared to the term. 
The mass mo of the 19th meson, the dilaton, remains finite above the phase 
transition so that Pmu —^ 0, and 

TT^ 

Pone loopi^dilcttOTl^ ^ P 

yu 


where only one degree of freedom for the gluon is included. 

A second step in the description of the gluon thermodynamics would consist 
in including a temperature-dependent gluon potential V^(T) as stipulated in 
[22,24], which could yield the right asymptotic behaviour (32). 

Fig. 7b is the equivalent of Fig. 7a for /i„ = = 250 MeV. These values of 

the chemical potentials correspond to a strong coupling between (uu) and (ss) 
with the set (11) while they correspond to a weak coupling with the set (12). 
Even at high temperature, it is no more possible to perform a linear fit of 
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the results. Instead, the chemical potentials introduce a dependence in the 
pressure. In the chiral limit (m„ = = 0), one would have 


P 


7 n - Nc 0^0 NcU"^ 

— 


60 


B 


(36) 


with ji = Hu = fig- Due to ^ 0, Eq. (36) has to be modified. However, the 
prescription used in the case /r = 0 (see Appendix A) does not work here if 
Hi > rui. From our results, one can see that the term is not modified by 
msi as expected, and the term is slightly decreased. 

The temperature behaviour of the energy density (21), (24) is exhibited in 
Fig. 8a, for Hu = l^s = 0. At high T, one can write 


7 


£ideal gas ~ —NcTT^T"^ - + 


20 -- - - * 


Nc 4 , rrig 
m. In —- 
ttT 


+ 


Nc 






IfivT^ 

while the obtained results seem to show a linear behaviour 

£ 10.269 + 1.620 10■^ 

£ 10.271 + 0.877 10“^ 


(37) 


(38) 

(39) 


for sets (11) and (12), respectively. When a fit is performed with the terms 
written in Eq. (37) plus an adjustable constant C, one obtains ~ 208 
MeV for set (11) and ~ 182 MeV for set (12), in agreement with the 
results obtained from the T behaviour of the pressure. 

The first order phase transition and the strong coupling between {uu) and 
(ss) (in the case of the set (11)) is well visualized by the sudden change of 
the energy value at Eg. For set (12), the weak coupling between the quark 
condensates and the presence of a second order phase transition for the u- 
quark (the s-quark still feeling a first order transition) is reflected by a slope 
breaking in e{T) at T„c and a small jump at Tgc- 

Fig. 8b is equivalent of Fig. 8a in the case of a dense strange matter Hu = 
Hs = 250 MeV. Here again, the chemical potentials introduce a dependence 
which writes with /r = : 

£ = ^NcTt^T^ + ^NcH^T^ + + B (40) 

zO Z 4 TT^ 
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in the chiral limit. From our results, one sees that rUs ^ 0 only slightly affects 
the second term of (40). It also introduces an additional constant term. 


For a pure ideal gas, the entropy density behaves as at /r = 0. When the 
non-vanishing s-quark mass is taken into account, one can use Eqs. (22), (27) 
and (37) to get the three first terms of the expansion of Ts in power of T^. 
One then has for s : 


7 Nr O 

S ^—N^T^ - ^mlT + 

15 6 


Nc 


4r7-i —1 


+ 


(41) 


where the dots represent more negative powers of T. However, the results of 
Fig. 9a show that one can approximate (41) by a linear expression : 

s 13.868 - 0.136 10"\ (42) 

s 13.839 - 0.912 10"^ (43) 


for sets (11) and (12) respectively. The coefficients of are not exactly equal 
to what would be obtained from combining Eqs. (29), (30), (38) and (39), 
although the numerical results are perfectly verifying Eq. (19). Similarly, the 
constants in Eqs. (42) and (43) are departing numerically from the second 
and third terms of Eq. (41). This reflects the limited meaning of simple fits 
like Eqs. (29), (38) and (42) in a restricted range of T^. This remark would 
apply to any calculation. It would then be dangerous to try to extract a bag 
constant from for instance P versus curve by a simple fit. 

Eig. 9b shows the temperature behaviour of s in the case of non-vanishing 
chemical potentials. Using (22), (36) and (40) together with the density after 
chiral restoration 

PH = ^ (44) 

one gets for an ideal gas : 

Ts = (45) 

expression which should be modified in order to take 7^ 0 into account. 
We observed that this non-vanishing only slightly affects the value of the 
factor multiplying T^. It should be noticed that above the transition, the 
entropy shows that all the available degrees of freedom are excited. However, 
as noted in Refs. [11,26], the massless Stefan-Boltzmann limit of QCD is not 
reached, because not all the gluon degrees of freedom are manifest in this 
model, and because the mass of the strange quark is not yet negligible. 
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5 CONCLUSION 


We have investigated the properties of the chiral transition in an infinite sys¬ 
tem of quarks with the help of a SU(3) flavour NJL model, built to account for 
the axial and scale anomalies. There is some freedom in chosing the parameters 
of the model. We illustrated our calculations for two choices, corresponding to 
large and moderate constituent quark masses. With the first choice, the tran¬ 
sition is always of first order. With the second choice, the transition can be of 
second order for small u and s chemical potentials. This is to be contrasted 
with many other NJL models (see e.g. Refs. [16,18]), where the precise form 
of the Lagrangian, the way of regularization and the choice of the parame¬ 
ters hardly allow for a first order transition. Let us mention that the chiral 
transition for the flavour SU(3) case is also studied in ref. [28] in a cr-model 
and in ref. [29] in a NJL model. These works however concentrate on the zero 
chemical potential case and on the modification of the meson masses. They 
include ^ [28] and ^ [29] terms. 

An important feature of our results lies in the possibility of having coupled, i.e. 
simultaneous, chiral transitions for light and strange quarks. These transitions 
tend to become uncoupled if either the chemical potential /r„ or the chemical 
potential is large enough. We remind that, in our model, the u and s 
condensates are only coupled indirectly through their coupling to the gluon 
condensate. As a consequence the latter is always modified in a quark (u or 
s) chiral transition. 


The scaled NJL model studied here allows the possibility of having first order 
transitions, with important jumps of the thermodynamical quantities. This 
is an interesting feature in view of the search of the quark-gluon plasma in 
heavy ion collisions. However, as we explained in section 3, the model should be 
made more realistic through, e.g. the introduction of vector mesons. It should 
be checked how the discontinuities will be modified by this introduction. 


We also investigated the equation of state for the chiral phase. The latter 
is very close to the relativistic free gas one. Departures from the massless 
Stefan-Boltzmann limit comes from the non-vanishing mass of the s quark 
in the chiral phase and from the “bag” constant, which can be identified as 
due to the chiral symmetry breaking of the vacuum (regardless of the non¬ 
zero current mass of the strange quark). Next steps in our investigations would 
involve the inclusion of the mesons degrees of freedom, which can be viewed as 
the quantum fluctuations of the quark fields [17], the study of their influence on 
the equation of state, especially in chiral asymmetric phase and the possibility 
of including all the gluon degrees of freedom. 
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Appendix 


A High temperature expansions for a hot ideal gas 


We give explicitly the expansion of thermodynamical quantities at zero chem¬ 
ical potential and high temperature. 

a) Pressure 


At high T, the pressure of an ideal gas of massless u quarks and massive s 
quarks is given by 


Pideal gas Pideal gas (^) T Pideal gas ('5) 


(AT) 


or 


Pi 


ideal gas 


= 2 - 




f 1, f fc" 

J ^ J ~E~i 


1 




■dk 


(A.2) 


with Eg = -|- rUg, rris being the bare s-quark mass. The first term of the 

r.h.s. of (A.2) has the usual behaviour 


Pideal gas(u) = 


7 


(A.3) 


The second term can be estimated in the following way : 

2]V. 

Pideal ^ hm / - —dv 

laeai gas\ ) ^^2 ^ ^^0 J I ^ 

1 


(A.4) 


where the infinitively small positive quantity e allows to regularize the T in¬ 
dependent term of Pideal gas{s) (sce below). One then has : 
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ideal gas 




n+1 ^—ne 


n=l 

oo 


oo 

/(.-I) 


3/2 


—nrus 




lim V {'^rrisd) 


TT^ 


(A.5) 


n=l 


Writing the modified Bessel function K 2 (z) as an ascending series in power 
of 2 [27] 


K 2 {z) = 2z 


-2 


1 - 


1 2 (fc + 1) + (A; + 3) / z"^ 

® k\{2 + k)\ 


k=0 


In IJ2 {z) 


(A.6) 


with ijj being the Digamma function and with the other modified Bessel func¬ 
tion h^z) 




(A.7) 


one may write (A.5) as a double sum which is absolutely convergent, owing 
to the convergence factor One can easily control the first terms of the 

expansion 


Pideal gas 


(s) = P4T^ -\-P 2 T‘^ + Pr In T 

+ Po +Pi- 2 T~'^\n T + p_ 2 T~^ 


(A.8) 


One has 




00 „-r. 


e-n£('_l)n+l ^ 


P4 = 


7T^ e^O 


lim y 


n+1 




E 


n=l 

00 ^-nef ^\n 


7r2 


= E_/Vc7r^ (A.9) 
180 ^ ^ 


Nc 2r ^e--(-l) 
P 2 = hm y-^— 


1 N 

= = (A-io) 


V 2 ~ 


1 rp -V 4 n. / \n —ne 1 'IZUlgP 

Pi in T = — Kmz hm > ( —) e In - 

47r2 ^ 2 


Nc 4 

= 

47r^ 

Nc 4 
47r^ ^ 


m /? °° °° 

In ^ lim y + hm y In n 
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n=l 
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.1 In + hm y y 
2 2 £^o da; 


a;=0 
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(A.ll) 


Nc 4 
= -r^^s 

Nc 4 
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4tt^ ^ 


1 rUsP d p 

— In-— 

2 2 dx 


(l_2i+-)c(-x) 


X=0 


1 rrisfd , 1 , tt 

in -H - m — 

2 2 2 2 


Nc 4, ms(3 
tTU, in-, 


^2 


StT ^ TT 
where ( is the Riemann ( function, 


Po = ^rn^ime^J2(-'^T (1) + ^ (3)} e" 


n=0 


1 


Nc 4/ „ I r 

= —27 + - lim 

Stt^ H 2A-ol + e-^ 
where 7 is the Euler constant, 
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^mlT-\\-2), 


487r" 


(A.13) 


as the Riemann function C(s) vanishes for s = —2. The p_2 term vanishes, 
since it involves the sum 


= C (-2) = 0. 

n=l 


(A.14) 


With (A.3), (A.8)-(A.14), the high temperature behaviour of the ideal gas 
pressure writes 
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ideal gas 


= —NcTt^T^ - —mlT'^ - 
60 12 " 


Nc 4 , 

— kTYI^ in — 
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+ 22±. 

ttT IGtt^ 


— 27) m 


7N, 
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(A.15) 


b) Energy density 
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Using the same technical procedure as before, one finds : 


^ideal gas ^ideal gas (^) “I” ^ideal gas 




(A.16) 
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With [27] 

A'l (j) 


~z~4h, WWi 1^1 


(A.18) 


(A. 19) 


and 


2,§,H(fc+l)!U 


2 \ ^ 



(A.20) 


one finds : 


^ideal gas (s) — + £2^"^ + £^ In T + £9 

+ £^_2T“Mn T + £_2r“^ + ... 


with 


7 

£4 = 3p4 = — NcTt '^ 
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Nc 2 


(A.21) 
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(A.23) 
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(A.24) 


In T = 3p(^ In T + lim ^ 
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where we have used (A. 11), 
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The term £_2 vanishes because it is proportional to C, (—2). 
results, one obtains : 


^ideal gas 
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c) Entropy density 


Using (22), (A. 15) and (A.27), one gets : 


T* = Um f; (’'“•'i* 

n=l 


nmsP 


s = —N^T^ - —mlT + ^mlT-^ + — (- 

15 G * Stt^ " 24 7r2 * ^ ^ 


(A.25) 


(A.2G) 
Gathering the 
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(A.28) 
2) + ...(A.29) 
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Figure captions 


Figure 1. 

Chiral phase transition Hsc versus /r„c for T = 0 MeV (a), T = 100 MeV (b), 
T = 161 MeV (c). The free parameters of the model are = 600 MeV, 
Xo = 125 MeV. The line with the crosses (open circles) corresponds to a 
vanishing ss [uu] condensate. 

Figure 2. 

Temperature behaviour of the constituent u- and s-quark masses and of the 
gluon condensate for = /i^ = 250 MeV (a) and = 200 MeV, fJ^s = 0 MeV 
(b) with = 600 MeV, Xo = 125 MeV. 

Figure 3. 

Chiral first order phase transition (Tc, fiuc, fJ^sc) ■ artistic view. 

Figure 4. 

Temperature behaviour of the constituent u- and s-quark masses and of the 
gluon condensate for Hu = IJ^s = ^ MeV (a) and Hu = = 250 MeV (b) with 

MO = 300 MeV, Xo = 80 MeV. 

Figure 5. 

Chiral phase transition fisc versus Hue for T = 0 MeV (a), T = 100 MeV (b), 
T = 130 MeV (c) and T = 165 MeV (d), with = 300 MeV, Xo = 80 MeV. 
Same conventions as for Fig. 1. 

Figure 6. 

Chiral phase transition {T^ Hue, Hsc), summarizing the results presented in 
Fig. 5 : artistic view. 

Figure 7. 

Equation of state for Hu = Hs = ^ MeV (a) and Hu = Hs = 250 MeV (b). 
Approximated linear fits above Tc are indicated in (a) while quadratic fits 
in are required for part (b). Results are shown for = 600 MeV, Xo = 125 
MeV (crosses) and M° = 300 MeV, Xo = 80 MeV (full dots). 

Figure 8. 

Energy density versus for Hu = Hs = ^ MeV (a) and Hu = Hs = 250 MeV 
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(b). Same conventions as in Fig. 7. 

Figure 9. 

Entropy density versus for = 0 MeV (a) and Hu = IJ^s = 250 MeV 

(b). Approximated linear fits, above Tg are indicated in (a) while quadratic 
fits for Ts are required for part (b). Same conventions as in Fig. 7. 
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